o ' Stability of undercompressive viscous shock 

profiles of hyperbolic-parabolic systems. 



< 



MOHAMMADREZA RA00FI*AND KeVIN ZuMBRUN^ 

February 2, 2008 

Abstract 



Extending to systems of hyperbolic-parabolic conservation laws 
results of Howard and Zumbrun for strictly parabolic systems, we 
show for viscous shock profiles of arbitrary amplitude and type that 
. necessary spectral (Evans function) conditions for linearized stability 

I established by Mascia and Zumbrun are also sufficient for linearized 

' and nonlinear phase-asymptotic stability, yielding detailed pointwise 

. estimates and sharp rates of convergence in L^, 1 < p < oo. 

o 

^ ' 1 Introduction 

g ■ Consider a (possibly) large-amplitude viscous shock profile, or traveling-wave 

• • ■ solution 
> 

(1.1) u{x — st); lim u{x) = u±, 

of a system of partially or fully parabolic conservation laws 

(1.2) ut + F{u)^ = {B{u)u^)^, 
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X G R, M, F e M", 5 G M"''", where 

(-) ^-if-)^ 

positive integer, possibly n (full regularization), 

and 

Re ^(62) > ^ > 0. 

Here and elsewhere, a denotes spectrum of a matrix or other linear operator. 
Working in a coordinate system moving along with the shock, we may without 
loss of generality consider a standing profile s = 0; we take s = from 
now on. 

Following [Z2], we assume that, by some invertible change of coordinates 
u — > w{u), followed if necessary by multiplication on the left by a nonsingular 
matrix function S{w), equations (11.21) may be written in the quasilinear, 
partially symmetric hyperbolic-parabolic form 

(1.4) A^wt + Aw, = (Bw,), + G, w= (^^'jj^ , 

G M""'^, w^^ G W, X G M, t G IR+, where, defining w± := w{u±): 
(Al) A{w±), All, A^ are symmetric, > 0. 

(A2) Dissipativity: no eigenvector of dF{u±) lies in the kernel of B{u±). 
( Equivalent ly, no eigenvector of A{A^)~^{w±) lies in the kernel of B{A^)~^{w±).) 
(A3) 5 = ( ° H , G = ( ° ) , with Reb{w) > 9 for some > 0, for all w. 



bj' 

and g{w,,Wo,) = 0{\w,\'^). 

Here, the coefficients of (11.41) may be expressed in terms of the origi- 
nal equation (11.21) . the coordinate change u — > w{u), and the approximate 
symmetrizer S{w), as 

io := S{w){du/dw), A := S{w)dF{u{w)){du/dw), 
^^'^^ B := S{w)B{u{w)){du/dw), G = -{dSw,)B{u{w)){du/dw)w,. 

Alternatively, we assume, simply, 

(Bl) Strict parabolicity: n = r, or, equivalently, ^a{B) > 0. 
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Along with the above structural assumptions, we make the technical hy- 



(HO) F, B,w, S e C^. 

(HI) The eigenvalues of := Aii{A^j)~^ are (i) distinct from 0; (ii) of 

common sign; and (iii) of constant multiplicity with respect to u. 
(H2) The eigenvalues of dF{u±) are real, distinct, and nonzero. 

(H3) Nearby -u, the set of all solutions of fll.ll) - fll.2p connecting the same 

values u± forms a smooth manifold {u*^}, S eU G M^, = u. 

Remark 1.1. Structural assumptions (Al)-(A3) [alt. (Bl)] and techni- 
cal hypotheses (H0)-(H2) admit such physical systems as the compressible 
Navier-Stokes equations, the equations of magnetohydrodymics, and Slem- 
rod's model for van der Waal gas dynamics \L2\ IZ3j . Moreover, existence of 
waves u satisfying (H3) has been established in each of these cases. 

Definition 1.2. An ideal shock 



is classified as undercompressive, Lax, or overcompressive type according as 
i — n is less than, equal to, or greater than 1, where i, denoting the sum of 
the dimensions i- and z+ of the center-unstable subspace of dF{u-) and the 
center-stable subspace of dF{u^) , represents the total number of character- 
istics incoming to the shock. 

A viscous profile f ll.ip is classified as pure undercompressive type if the 
associated ideal shock is undercompressive and i = 1, pure Lax type if the 
corresponding ideal shock is Lax type and i = i — n, and pure overcompressive 
type if the corresponding ideal shock is overcompressive and i = i — n, i as 
in (H3). Otherwise it is classified as mixed under-overcompressive type; see 



Pure Lax type profiles are the most common type, and the only type 
arising in standard gas dynamics, while pure over- and undercompressive 
type profiles arise in magnetohydrodynamics (MHD) and phase-transitional 
models. Mixed under-overcompressive profiles are also possible, as described 
in |LZ2l IZHj . though we do not know a physical example. In the pure Lax 



potheses: 



(1-6) 




M_ X < St, 

X > st, 
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or undercompressive case, {u^} = {u{- — S)} is just the set of all translates 
of the base profile u, whereas in other cases it involves also deformations of 
u. For further discussion of existence, structure, and classification of viscous 
profiles, see, e.g., [LZ2l [Ml [MZ2l lMZ3l [MM [ZH [Z2l [Z3], and references 
therein. 

Definition 1.3. The profile u is said to he nonlinearly orbitally stable if 
u{-, t) approaches u^^*^ as t ^ oo, as defined in (H3), for any solution u of 
M.2\) with initial data sufficiently close in some norm to the original profile 
u. If, also, the phase 6{t) converges to a limiting value 5*, the profile is said 
to be nonlinearly phase- asymptotically orbitally stable. 

An important result of |MZ3j was the identification of the following stabil- 
ity criterion equivalent to linearized orbital stability of the profile, 
p > 1, where -D(A) as described in |GZt IZHj denotes the Evans function asso- 
ciated with the linearized operator L about the profile: an analytic function 
analogous to the characteristic polynomial of a finite-dimensional operator, 
whose zeroes away from the essential spectrum agree in location and multi- 
plicity with the eigenvalues of L. 

(P) There exist precisely i zeroes of D{-) in the nonstable half-plane 
MA > 0, necessarily at the origin A = 0. 

As discussed, e.g., in |MZ4t [Z2| IZSj . under assumptions (Al)-(A3) and 
(H0)-(H3), (V) is equivalent to (i) strong spectral stability, cr(L) C {MA < 
0} U {0}, (ii) hyperbolic stability of the associated ideal shock, and (iii) 
transversality of m as a solution of the connection problem in the associ- 
ated traveling-wave ODE, where hyperbolic stability is defined for Lax and 
undercompressive shocks by the Lopatinski condition of |Mll IM2^ IM3I, IFrej 
and for overcompressive shocks by an analogous long-wave stability condition 
[Zlt IZ2] . Here and elsewhere a denotes spectrum of a linearized operator or 
matrix. 

The stability condition holds always for small-amplitude Lax profiles 
[HuZt IPZt IFreSj . but may fail for large-amplitude, or nonclassical over- or 
undercompressive profiles |AMPZlt \GZ\ \FreZ\ \ZS\ [Z2] . It may be readily 
checked numerically, as described, e.g., in [Brit IBr2l IBr3l [BrZt IBDG] . It was 
shown by various techniques in [MZTl [MZ2l [MZ3l MM, MZ5l [Ral [HRl [HRZ] 
that the linearized stability condition {T>) is also sufficient for nonlinear or- 
bital stability of Lax or overcompressive profiles of arbitrary amplitude. In 
the strictly parabolic case (Bl), this result was extended in [HZ] to shocks 
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of arbitrary amplitude and type, in particular to shocks of under- or mixed 
over-undercompressive type. However, up to now, it had not been verified 
for under- or under-overcompressive profiles of systems with real viscosity. 

In this paper, we establish for shocks of any type and general systems 
satisfying (Al)-(A3) [alt. (Bl)] and (H0)-(H3) that {V) is sufficient for 
nonlinear phase-asymptotic orbital stability. More precisely, denoting by 

(1.7) af < af < ■ ■■ < 

the eigenvalues of the limiting convection matrices A± := df{u±), define 



:i.8) 9{x, t) := ^ (1 + t)~^/\-\--Jt\VLt + ^ (1 + t)-V2g- 



\x-aJt\^/Lt 

5 



a7<0 aj>0 



(1.9) 

and 
(1.10) 



ipi{x,t) ■.= x{x,t) ^(1 + |x| +t)-^/2(l + |x-a7t|)-i/2 

a-<0 

+ x(x, t) ^ (1 + + t)-^/\l + \x- ap 



1/2 



aT>0 



il^2{x, t) := (1 - xix, t))(l + \x- a^tl + ti/2)-3/2 
+ (1 - X(x, t))(l + |x - a+t| + tV2)-3/2^ 



where xi^y t) = 1 for x G [a^^t, a^t] and zero otherwise, and L > is a 
sufficiently large constant. 

Then, we have the following main theorem. 

Theorem 1.4. Under assumptions (A1)-(A3) [alt. (Bl)], (H0)-(H3), and 
{V), the profile u is nonlinearly phase- asymptotically orbitally stable with 
respect to H'^ initial perturbations Uq, with ||(1 + |xp)~^/^Mo(x)||j|^5 < Eq 
sufficiently small. More precisely, there exist 6{-) and 5* such that 

\u{x, t) - (x) I < CE^{Q + ^1 + ^/-s) (x, t), 

\d,{u[x,t)-u''^''^'\x))\ < CEoie + A+^2){x,t), 
(1.11) < CEo, 

m\<CEoil + tr\ 
m\<CEo{l + t)-'/\ 
where u denotes the solution of fll.2p with initial data Uq = u + Uq. 
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In particular, Theorem 11.41 yields the desired result of nonlinear stabil- 
ity in the undercompressive or mixed case, effectively completing the one- 
dimensional stability analysis initiated in [ZHt IMZ3] . 

Remark 1.5. Pointwise bound fll.lip yields as a corollary the sharp U' decay 
rate 

(1.12) \u{x,t) -u^'^^^'^{x)\lv <CEo{l + tY'^^^-l\ l<p<oo. 

The main new difficulty in the analysis beyond those faced in the strictly 
parabolic case of |HZ| is to control higher derivatives in the absence of 
parabolic smoothing. We accomplish this by a modification of the fixed- 
point iteration scheme introduced in [HZ] , changing to an implicit iteration 
scheme in order to avoid loss of derivatives as discussed in Remark 16.21 

This is a standard device in situations of limited regularity, especially for 
quasilinear hyperbolic equations. However, here the situation is complicated 
by the nonlocal character of the defining integral equations, which appears to 
prevent the standard treatment of regularity by energy estimates. The simple 
resolution to this problem is that, by appropriate choice of implicit scheme, 
we obtain a nonlocal system of integral equations that admits also a local 
description in terms of a symmetric hyperbolic-parablic system amenable to 
the same type of energy estimates used to control regularity in the study of 
Lax and overcompressive shocks in [MM [Z2l[Ra] . 

New physical applications beyond those of |HZ] are to undercompressive 
waves in MHD and, with slight modification following |LRTZ] . to weak det- 
onation waves in reactive compressible Navier-Stokes equations. The latter 
we intend to treat in a future work. 

Plan of the paper. In Sections |2] and [3l we recall the basic profile 
bounds and linearized estimates obtained in |MZ3l IZHl IHRZ] , and in Section 
m the convolution estimates established in [HZj . In Section [5], we recall (a 
slight modification of) an auxiliary energy, or "hyperbolic-parabolic damp- 
ing" estimate, established in |MZ4^ IZ2t IRaj . along with a more standard 
weighted estimate. In Sections [6l [71 and [8], we carry out the main work of 
the paper, introducing a crucial implicit version of the iteration scheme de- 
scribed in |HZj . establishing local existence by energy estimate, and then 
showing by a combination of estimates like those of [HZj and [Z2j that this 
is contractive in an appropriate norm encoding the claimed rates of decay. 
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2 Profile facts 



We first recall the profile analysis carried out in |MZ3t IHRZj , a slight gener- 
alization of Corollary 1.2, [ZHJ, which in turn generalizes results of [MP] in 
the strictly parabolic case. Profile u{x) satisfies the standing-wave ordinary 
differential equation (ODE) 

(2.1) B{u)u' = F{u) - F{u^). 
Considering the block structure of B, this can be written as: 

(2.2) F\u',u'') = F\ui,u'J) 
and 

(2.3) hiu'Y + b2{u''y = F'\u\ u'') - F'\ui, u'l). 

Lemma 2.1 f pvIZ8l IHRZ] ). Given (H1)-(H3), the endstates u± are hyper- 
bolic rest points of the ODE determined by l{2.3\) on the r -dimensional man- 
ifold h2.2^) . i.e., the coefficients of the linearized equations about u±, written 
in local coordinates, have no center subspace. In particular, under regularity 
(HO), 

(2.4) DiD\{u\x) - u±) = 0(e-"l^l), a > 0, < j < 5, z = 0, 1, 
as X ±oo. 

Proof. By (HI), (12. 2p may be solved for = h{u^^), reducing the problem 
to an ODE in u^^. Under assumptions (H1)-(H3), this is a nondegenerate 
ODE of which u^J are hyperbolic rest points; see |MZ3i[Z2] . The family is 
thus the intersection of the unstable manifold at m_ with the stable manifold 
at u+, both of which are by (HO) and standard invariant manifold theory. 
This intersection is transversal as a consequence of (T>), |MZ3] . hence is 
in 6 by the Implicit Function Theorem, and in x by (HO) and the 
defining ODE. Finally, (12. 4p follows from hyperbolicity of u±, by standard 
ODE estimates on the ODE and its variations about u^. For further details, 
see jMZ3l [Z2]. □ 
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3 Linearized estimates 

We next recall some linear theory from |MZ3t IZHj . Linearizing (11.21) about 



u 



■), (5* to be determined later, gives 



(3.1) 
with 



Vt 



(3.2) B^*{x) := B{u^'{x)), A^*{x)v := dF{u^'{x))v - dB{u^' {x))vui* . 

i?(±oo), and considering Lemma l2.ll it 



Denoting := A{±oo), B^ 
follows that 



(3.3) 



A^'(x)-A- 



0(e 



'TI\X\ 



\B^'{x)-B- \ = 0(e 



-ri\x\ 



as x —oo, for some positive r). Similarly for and B~^, as a; — > +oo. Also 
1^4''* (x) — A^l and \B^*{x) — are bounded for all x. 

Define the (scalar) characteristic speeds < ■ ■ ■ < {as above) to be 
the eigenvalues of A^, and the left and right (scalar) characteristic modes if, 
rf to be corresponding left and right eigenvectors, respectively (i.e., A^rf = 



etc.), normalized so that If ■ 



61 and Ij ■ 



61. Following 



Kawashima [Kawj . define associated effective scalar diffusion rates (3f : j 
1, ■ ■ ■ , n by relation 



(3.4) 



0\ 







diag L^B^R^, 



where := {if, . . . , /^)*, := {rf, . . . , rf) diagonalize A^. 
Assume for A and B the block structures: 



A 



All Ai2 
A21 A22 



B 





B21 B22 



Also, let a*{x), j = 1, . . . ,{n — r) denote the eigenvalues of 



A^ := All — ^12-822^-821, 
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with r*(x) G ^ associated left and right eigenvectors, normahzed so 

that Ifvj = 61- More generally, for an m*-fold eigenvalue, we choose (n — r) x 
m* blocks L* and R* of eigenvectors satisfying the dynamical normalization 

Lfd^R* = 0, 

along with the usual static normalization L*^Rj = ^ilm*] as shown in Lemma 
4.9, |MZlj . this may always be achieved with bounded L*, R*. Associated 
with L*, R*, define extended, n x m* blocks 



Eigenvalues a* and eigenmodes C*, 71* correspond, respectively, to short-time 
hyperbolic characteristic speeds and modes of propagation for the reduced, 
hyperbohc part of degenerate system fll.2p . 

Define local, mj x rrij dissipation coefficients 

r]*{x) := -LfD,R*ix), j = 1, . . . , J < n - r, 

where 



Djx) := AuB^i 



A21 — A22B22B21 + A^B22B2l + -822 -822^-^21) 



is an effective dissipation analogous to the effective diffusion predicted by 
formal, Chapman-Enskog expansion in the (dual) relaxation case. 

The Green distribution (fundamental solution) associated with (13.11) is 
defined by 



(3.5) G{x,t;y) := e^'''5y{x). 

or, equivalently, 

Gt~L^*G = 0, \imG{x,t;y)=6y{x). 

Recalling the standard notation errfn(z) := ^ e~^^d^, we have the fol- 
lowing pointwise description. 
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Proposition 3.1 ( [MZ3j ). Under the assumptions of Theorem I.4, the Green 



distribution G{x, t; y) associated with the linearized equations fl3.1l) may he 
decomposed asG = H + E + G, where, for y < 0: 



J 



H{x,t;y) := a*-l(x)a*(y)7^*(x)C;(z/, (l/) 
(3.6) 7' 

= ^7^*(x)0(e-^«*)^._,.,(-^/)£**(2/), 



where the averaged convection rates a* = a*{x,t) in liS. 6\) denote the time- 
averages over [0,t] of a*{x) along backward characteristic paths z* = z*{x,t) 
defined by 

dz* /dt = a*{z*), z*{t) = X, 

the dissipation matrix (* = (*{x,t) G ]R'"^^™'j is defined by the dissipative 
flow 

dc;/dt = -r]*{z*)q, c;(o) = /„.,, 

and 6x-a*t denotes Dirac distribution centered at x — a*t. 

(3.7) ^(^'^^5^/) = ^^^ ^i(2/,t), 

j=i \s=d 



for y < and symmetrically for y > 0, with 
(3.9) \ltk\<C^ \{d/dy)l%\<Ge-^\y\, 
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and 
(3.10) 

|a:,,G(x,t;y)|<Ce-''(l^-^l+*) 



n 
k=l 



/2 ^-(x-y-altf /Mt ^-^x+ 



'J {t-\yhi\)?/Mt^-r,x+ ^ 



I V .-1/2 -{x-a 

a^>0,aj<0 



Oj. >0, aj>0 



< |a| < 2 for y < and symmetrically for y > 0, for some rj, C , M > 0, 
where are as in Theorem \1.4l > 0, denotes the positive/negative 
part of X, indicator function X{\a-t\>\y\} ^-^ /^'^ \^k^\ — Ivl ^ otherwise. 
Moreover, all estimates are uniform in the supressed parameter d. 

Proof. This is a restatement of the bounds estabhshed in |MZ3j for pure un- 
der compressive, Lax, or overcompressive type profiles; the same argument 
apphes also in the mixed under-overcompressive case. Also, though it was 
not explicitly stated, uniformity with respect to is a straightforward con- 
sequence of the argument. □ 



Corollary 3.2. Under the assumptions of Theorem I.4 and the notation of 
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Proposition \3.1\ 



\ej{y,t) — ej{y, +c>o)\ < Cerrfn( ^^^^ ^ ), some a > 



M^ft 



a->0 



+ Ce-^'l^l errjfi 



y + a^t 



errfn 



\dye,{y,t)-dye,{y,+oo)\ < Cr'/' 



|j/+a-t|VMt 



/or y < 0, anc? symmetrically for y > 0. 

Proof. Straightforward calculation using fl3.8p and (13.91) : see |MZ3j . 
From now on, let 



(3.12) 



/ei\ 



□ 



Corollary 3.3 (|HZj). Under the assumptions of Theorem \1.4\ and the no- 
tation of Proposition \3.1l 

/oo 
e{y,+oo))idu'/d5)isMdy = Ii 
-oo 

Proof. This follows from the standard fact that L^* {du^ j d5)\s, = 0, hence 
G{x,t; y) {du'/d5) {y) dy = idu'/d5) (x) 
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which, together with the fact that E = {du^ /dS)\s,{x)e{y,t)) represents the 
only nondecaying part of G{x,t;y) under stabihty criterion (V), yields 

r+oo 

{du'/d6)\sM / e{x,+oo;y){du'/d6)\sMdy = {du'/d6)\sM 



in the limit as t 



— oo 
-OO. 



□ 



Proposition 3.4 (Parameter-dependent bounds). Under the assumptions of 
Theorem \1.4\ and the notation of Proposition \3.1l 



(3.14) 



ds^H{x,t-y)=0{H) + dyO{tH), 



(3.15) 





1^5. e| 


< C\e\ 




\9s,et\ 


<C\et\, 


l^5,(e(i/,t) - 


e{y,+oo))\ 


<C|(e(y,t)-e(y,+oo))| 




\95,ey\ 


— ^I^J/ 1 5 




Cyiy, +oo) 


< C\{ey{y,t) - Cyiy, +oo] 




\d5,eyt\ 


< C\eyt\, 



and 
(3.16) 



\ds,d:^yG{x,t;y)\<C\d:^^ G{x,t;y)l 



< I a I < 2 for y < and symmetrically for y > 0, for some C > 0. 

Proof. These follow by the same argument used to establish the parameter- 
dependent bounds of Proposition 3.11, |TZ1] . using the additional fact that 
neither speed s nor endstates u± depend on the choice of 5* to obtain better 
decay estimates on certain terms. 

Bounds (13.141) and (I3.15P follow by direct calculation, together with the 
observations (obtained similarly as bounds established in the proof of Propo- 
sition 3.11, [TZlj . using parameter-dependent asymptotic ODE bounds) that 



ds,a*, ds,a*, ds.TZ}, ds,Cj 



0(1) 



and 

(3.17) d5^d^,l,k^{y) = 0{e-^\y^), < \a\ < 2. 

Bounds (I3.16P follow by the argument of |TZ1] . but using the fact that a = 
a^, (3^ (since u±) do not depend on 6^,, hence "bad" factors tds^a = 0{t) do 
not appear, but only factors 0(1) or better. □ 
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Remark 3.5. The additional factor t in the righthand side of fl3.14l) may 
be absorbed in time-exponential decay of H, i.e., tH obeys the same decay 
bounds as H, but with slightly smaller time-exponential decay rate. 



4 Convolution lemmas 

We shall make use of the following technical lemmas proved in |HZt IHR] . 
Lemma 4.1 (Linear estimates I). Under the assumptions of Theorem I.4 



(4.1) 



00 

00 



00 
+00 



e 



\G{x, t; y)\{l + \y\r^/' dy < 0(6 + tfj, + ^2){x, t), 
\e,{y,t)\{l+\y\r^/'dy<C{l + tr'/\ 
iy,t)\{l+\y\)-'/'dy<C, 

/+00 
\e{y,t) - e{y, +oo)|(l + \y\)-'/' dy < C{1 + ^)-^/^ 
-00 

for < t < +00, some C > 0, where G and e are defined as in Proposition 

EH 

Lemma 4.2 (Nonlinear estimates I). Under the assumptions of Theorem 
(4.2) 

ft /•+00 

\Gy{x, t - s; y)\^{y, s) dyds < C{e + Vi + ^2)(a;, t), 



'0 J-OD 

"t-l r+oo 



/ / \G^y{x,t- s;y)\<i/{y,s) dyds <G{9 + '4^1 + '^2)ix,t), 

Jo J -00 

pt P+OO 

/ / \eytiy,t-s)\<il{y,s)dyds<G{l+t)-\ 

Jo J -00 

r+oo r+oo 

/ / \ey{y,+cx^)\m{y,s)dy<G{l+t)-^/\ 
Jt J -00 

ft r+oo 

/ / \ey{y,t-s)-ey{y,+(x^)\m{y,s)dyds<G{l+t)-^/'', 

Jo J -00 
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for < t < +00, some C > 0, where G and e are defined as in Proposition 
m\and 



^{y, s) := (1 + sf'\s-^/\e + ^1 + ^2f{y, s) 

+ (1 + S)-1(0 + ^1+V^2)(y,s). 



(4-3) 



Lemma 4.3 (Nonlinear estimates II). Under the assumptions of Theorem 



1.4 




(4.4) 



t f + OD 

\Gy{x, t - s; y)\^iiy, s) dyds < C{9 + ijji + V'2)(a;, t), 

J -00 

t— 1 P+CO 

/ \G^y{x,t- s;y)\^i{y,s)dyds < ^(0 + + V^2)(a:, t), 
J -00 

t P + OO 



/ / \eyt{y,t-s)\^^{y,s)dyds<C{l+t)-\ 

Jo J -00 

pt P+OO 

/ / \eyiy,t-s)\^,iy,s) dyds < Gil +t)-'^ 

Jo J -00 



and 
(4.5) 



a +00 
\G{x, t - s; y)\<l>2{y, s) dyds < G{9 + t/^i + ^2){x, t), 
-00 

pt p+00 

/ / G^^2{y,s)dyds<G{e + i)i + i)2){x,t), 

Jo J -00 

ft p+oo 

-3/2 




Pt p+00 

/ / \etiy,t-s)\^2iy,s)dyds<Gil + ty 

Jo J -00 

\e{y,t-s)-e{y,+oo)\<!>2{y,s)dyds < C(l + t)" 



J~oo 



for < t < +00, some G > 0, where G and e are defined as in Proposition 
I3l\and 

<l>i(y, s) := e-''\y\s-'/'i9 + + ^2)iy, s) < Ge-^^y^^^-'/^l + s)-\ 
<l>2{y,s) :=e-''l^l(l + s)~=^/2_ 
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Lemma 4.4 (Linear estimates II). Under the assumptions of TheoremVLj 

3 

if l'^^o(a;)|, |<9a;t'o(a;)| < Eq{1 + 1x1)^2, Eq > 0, then, for some 6' > 0, 



(4.7) 



oo 

+ 00 



H{x,t- y)vo{y)dy < CEoe-'\l + |a;|)-t 
H,{x,t;y)vo{y)dy < CEoe-'\l + \x\)--2 



□ 



and so both expressions are dominated by Eo{ipi + 4'2)- 

Proof See jHR |[HR2] . 

Lemma 4.5 (Nonlinear estimates III). Under the assumptions of Theorem 
if |T(y,s)| < s-^'^{:ipi + ^2 + d){y, s) + s-^/^e-'^\y\ then 



t r+oo 




H{x,t- s;y)T{y,s)dyds < C^ipi + V'2)(a;,t), 



J-oo 



pt P+OO 

/ / {H^, - Hy){x,t - s;y)T{y,s)dyds < Ciipi +ilj2){x,t), 

Jo J-oo 



t P + OO 




\Gx{x,t- s;y)\T{y, s) dyds < C{tpi +ilj2){x,t). 



t—l J — oo 



Proof. The proof of the first inequality is very much similar to that of the 
similar estimates proved in [HRJ . Here we prove only the part which contains 
the convolution of H against s~^^^il)i. A typical term in (14.81) coming from 
ipi would be dominated by a term of the form 

C e-*(*-^)x(a; - d*{t - s), s)s--^{l + s)^^ 

X (1 + |x — a*(t — s) — a~s\)~'^ds 



C 



now, usmg 



-»?o(t- 



1 



'\{x - a*{t - s),s)s 2(1 + s) ^ 



< 



l+|a+b| — l+|a| 



X (1 + |x — t — (a* — ttj )(t — s)|) -^ds; 
the above would be smaller than 



C 



x(l + |a; - a-t\)~i{\l + (a* - a^){t - s)\)^ds. 
16 



Notice that {a* — ) < C . Now we use half of rj to neutrahze t — s, and to 
get 

C(l + |x - art|)-5 ^ e-^(*-^)x(a; - a*{t - s), + s^Us 

On the one hand, this is obviously dominated by 

C{1 + \x-a-t\y^{i+t)-^. 

which is absorbed by a factor of ipi, if a; G [a|ft,a^t]. On the other hand, 
x{x — a*{t — s),s) in the integral means that the above expression would 
vanish if a; > Mt, for some fixed M. Therefore the above expression is 
always less than a factor of ipi + ip2- 

The second inequality follows by an identical proof, using — Hy ~ H. 
The proof of the third inequality is similar, using the fact that 

C\Gy{t -s)\< e~'^'-'\t - s)-ie-l"~«'l'/^(*~^) 

for t — 1 < s < t; see |MZ2] for similar calculations. □ 



5 Auxiliary energy estimate 

We shall require the following auxiliary energy estimate adapted essentially 
unchanged from jMZ4l[Z2l[Raj . Let 

(5.1) ut + f{u), - {B{u)u,)^ = {du/d5)\s^{x)^{t), 



(5.2) 



u 



u 



Lemma 5.1 ( |MZ4t [Z2l IRaj ). Under the hypotheses of Theorem \1.4\ let up e 

H^^ , and suppose that, /or <t<T, the suprema of\5\ and I7I and the W^''^ 
norm of the solution u = {u^,u^^Y of (15.11) . (15. 2p each remain bounded by a 
sufficiently small constant C > 0. Then, for all < t < T , 



(5.3) \u{t)\l, 



<Ce-''\u{0)\ls + C 



^-'^^'-\\u\l + \S\' + h\'){T)dr. 
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Proof. This follows exactly as in the 7 = case treated in [MZ4t 
observing that term 

{du/d6)\s^{xh{t) 

is of the same form as terms 

{du/d6)\sM-m 

already arising in the nonlinear perturbation equations in the former case. □ 

We require also the following much cruder estimate adapted from [HRj . 

Lemma 5.2 ([HRj). Under the hypotheses of Theorem I.4 let Eq := ||(1 + 
|xp)~^/'*Uo(a;)||j^5 < 00, and suppose that, for < t < T, the suprema of \5\ 
and I7I and the W^'°^ norm of the solution u = {u^,u^^y of (15 .ip . (15.21) each 
remain bounded by some constant C > 0. Then, for all < t < T , some 
M = M{C) > 0, 

(5.4) 11(1 + \x\'y'/Mx, t) His < Me''' (e, + + \^\'){t) dr) . 

Proof. This follows by standard Friedrichs symmetrizer estimates carried out 
in the weighted norm. (Recall, these plus several more complicated esti- 
mates are used in the proof of Lemma [5. 1[ ) □ 



Remark 5.3. An immediate consequence of Lemma 15. 2[ by Sobolev em- 
bedding and equation (15. ip . is that, if Eq := ||(1 + \x\^)~^^^uq{x)'^h^, II'^IIhs, 
|(5(-)| and |7(-)| are uniformly bounded on < t < T, then 

|(l+|xp)-3/4^(x,t)|, |(l+|xp)-3/4„,(x,t)| 

are uniformly bounded on < t < T as well. 



6 Fixed-point iteration scheme 

We now introduce the fixed-point iteration scheme by which we shall simulta- 
neously construct and estimate the solution of the perturbed shock problem. 
Our starting point, similarly as in [HZ|, is the observation that 

(6.1) u{x,t):=u{x,t)-u^*+^^*\x) 
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satisfies 



(6.2) ut -L^'u = Q^* (u), + 5{t) {du^/d5) is, + R^' (5, u), + 5^* (5, 6, 



where 
(6.3) 

Q'5*(m,m,) := (^f{u^*) + A{u^*)u-fiu^*+u)'^ + (^B{u^* +u)-B{u^*)y, 

= 0{\u\'^ + \u\\u^\), 
Q^'{u,u^)x = 0{\u\\u^\ + + |m||mxz|), 



R^* := \^A{u^'{x)) - A{u^*+^^'\x))ju = 0{e~'^^^^\6\\u\), 

(6.4) i?f* = 0(e-^N|5|K|), 
Rt = Oie-^\MM), 

equation and 

(6.5) S'^ := 5(^idu'/d5)^s,+5it) - {du'/d5)^s,) = 0(e-''l^l|5||5|) 

so long as |m| remains bounded, by Taylor's Theorem together with (12.41) . 
Accordingly, for given (5"~^(-), define u"' to be the solution of 



hoo 

M"(x,t) = / iH^'' + G--')ix,t;y)ur\y)dy 

J — oo 

pt P+OO 

+ / / iH--' + G--')ix,t-s;y) 

Jo J -oo 

X 5^*""'(5"-\5"-i)(?/,s)dyrfs 

ft f + OO 

G;-\x,t-s;y) 

'0 J-oo 

n-l 



(6.6) 




X {Q^*"' (m", O + i?-^*"" {6''-\u'')){y,s)dyds 



t f + OO 

H^~\x,t-s;y) 




'0 ^-oo 

n-l 
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where i?" ^, G" ^ are the parts of Green distribution ^ of hnearized equa- 
tion around u^*" , and 



n-l 



(6.7) u"^-^ ^ uo - u^*^' . 

Further, set 

/oo 
(e-i(^,t)-e"-i(|/,+oo)K-^(|/)dy 
-oo 

(e"-i(y,t-5)-e"-i(l/,+oo)) 



and 



(6.9) 



oo 

t p+oo 



./-oo 

X 5^*""'((5"-\(5"-^)(y,s)dyds 
(6.8) + / / ie;-\y,t-s)-e--\y,+oo)) 

Jo J-oo 



+ / / e^-\y.+^)S'''"-\5^-\5^-'){y,s)dyds, 

J t J — OO 

/' + 00 Z' + OO 

- / / e;-i(y,+oo) 

J-oo 

X (Q^-"-'(k",<) +i?^*"-'(5"-\K"))(y,s)ciyds 



/ + 00 
e"-^(y,+oo)<-^(y)dy 
'OO 

+ / / e''-\y,+^)S'*"-\5^-\5''-'){y,s)dyds 

Jo J-oo 

p-\-oc r-\-oc 

- / / e;-\y,+oo) 

Jo J-oo 

X {Q'*''-\u'^,u:) + R'*''-\S^-\u^)){y,s)dyds. 
Define an associated iteration map T, formally, by 
(6.10) (5\5:)^T(5^-\5r'). 
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Lemma 6.1. Under 06.61) - fl6. 91) . := + -u"^*" ^ satisfies 

(6.11) u7 + mn.-{Bmu:). = {s-{t)-6-'\t))^\, 



oo 

5"(t) = - / er\y,tK-\y)dy 



— oo 

t f+ca 




(6.12) Jo J-oo 

»t /•+00 



+ / / 6^1(2/, t-.) 
Jo J-oo 



with initial data «"(■, 0) = Mo + (m'^'" — m"^*" ), anc? therefore satisfies 

f ll.2l) wt/i initial data Uq if and only if 

(6.13) (5«,5^) = (5-1,5-1), 

z.e., (5", 5") zs a fixed point ofT, in which case also 5(0) = 5(+oo) = 0. 

Proof. Equation (I6.12p follows immeditately upon differentiation of (16. 8p . 

From (16.81) . we obtain, further, that 5"'(+cxo) = 0, and 

(6.14) 



hoo 

5"(t)-5"(0)=- / e--\y,t)ut\y)dy 

e-'\y,t - s)S''"-\6"-\6"~'){y,s)dyds 

e--\y,t-s) 
X {Q''"-\u^,u:) + R''"'\6''~\un)(y,s)dyds. 



oo 

t p+oo 




'0 J-oo 

-t f+OO 




'0 J-oo 

n-l 



Setting t = +00 in (16.141) . and comparing with (16.91) . we find, therefore, that 
5«(0) = 5"(+oo) = if and only if = 5^^"^ 
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From fl6.6p and fl6.8p we conclude that 



i-oo 

J — oo 

/•t r+oo 

+ G''-\x,t- s;y)S^*"'\S''-\S''-^)iy,s)dyds 

Jo J -oo 
pt P+OO 

+ / / G--\x,t-s;y) 



(6.15) 



oo 

n-1 



X 



+ (5-(t)-5"(0))^|,.-i 



and thus, by Duhamel's Principle, 
(6.16) 

7x"(t) - = Q^*""'(m", + i?''*""'(5"-\M"),. + ^^*""'(5"-\5"-^) 

— „ „ rn — 1 . 

Setting = m^ + m^*" \ we then obtain (16. lip by a straightforward 
calculation comparing with (16. 2p . with the claimed initial data 

(6.17) ^"(■, 0) = + = ^0 + (^^."-^+^"-^(0) - u^-'). 

Note that the righthand side is equal to uq if and only if 5"~^(0) = 0, or, in 
case 5" = 5"^^ (a fixed point), if 5"(0) = 0, or equivalently 5" = □ 

Remark 6.2. Other than a slight notational change 6 ^ + 6 made 
to simplify the exposition, the difference between this iteration scheme and 
the one used in |HZ] in the strictly parabolic case is that we have made 
it implicit in m", i.e., u appears everywhere on the righthand side of the 
integral equations with index n rather than n — 1. By this change we preserve 
regularity properties, as encoded in the nonlinear structure of equation (16. lip : 
see Lemma Em By contrast, the explicit version of [HZj is not associated with 
a (favorable) nonlinear equation, and so would lose derivatives, preventing 
the iteration scheme from closing. 

Remark 6.3. Note that (I6.1ip - (l6.12p form a closed system for {u'^,5"'), in 
the form of a true Cauchy problem; that is, the values of (m", 5") at time T 
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depend only on values for < t < T, and not on future times. By 06.61) . we 
have, evidently, 



n-l 



(6.18) u"(-,0) = u[J 

7 Local existence 

Lemma 7.1 {H^ local theory). Under the hypotheses of Theorem \1.4\ let 
El := \\u^{x)\W. + + l<^*""'l < oo. 

Then, for T = T{Ei) > sufficiently small and C = C{Ei,T) > suffi- 
ciently large, there exists on < t < T a unique solution 

(m",5") G L^{H\x)-t) X C\t) 

of fl6TT|) - (16121) ■ satisfying 

(7.1) Iklln^t), |^"(t)| <Ci?i. 

Proof. Short-time existence, uniqueness, and stability follow by (unweighted) 
energy estimates in similar to (15. 4p combined with more straightforward 
estimates on 5" carried out directly from integral equation (I6.12p . using a 
standard (bounded high norm, contractive low norm) contraction mapping 
argument like those described in |Z2[ [Z3j . We omit the details. □ 



Remark 7.2. A crucial point is that equations f l6.1ip -( l6.12p depend only 
on values of (m", 5") on the range t G [0, T]; see Remark [6.31 

8 Proof of the Main Theorem 

We are now ready to prove the main theorem. Define norms 

(8.1) \h\B, := \h{t){l + tfl\o.^,) + |/i(t)(l +t)Uo.(,), 
and 

(8.2) I^Ib, := \ f{e + + ^2)"Vi.-(x,t) 
and Banach spaces 

(8.3) Bi := {h ■ \h\B, < +oo}, B2 := {g : \g\B2 < +00}. 
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Lemma 8.1. Under the hypotheses of Theorem 1.41 

G L^{t,H^{x)) X 
satisfy f l6.11l) -( l6.12p on < t < T, and define 

(8.4) at):= sup {(\u-\ + \u:\)ie + ^,+i;2r\x,s) + \S-is){l + s)\). 



x,Q<s<t 

IfC{T), \\uq~^\\h5, and \S'^^^\bi are bounded by Cq > sufficiently small, then, 
for some e > 0, (i) the solution (u", 5""), and thus ( extends toO<t<T + e, 
and (a) C bounded and continuous onO<t<T + e. 

Proof. By f l6.18p and Lemma [STTl smallness of ((T), \6"'~^{s)\bi, and 



^0 IIh^ = II«"(-,o)|| 



(recall ( 16.18P ) together imply boundedness (and smallness, though we don't 
need this) of ||m||h5 and \6^\l°° on < t < T. By Lemma 17. this implies 
existence, boundedness of ||m||j^5, ||5||vi/i,oo onO<t<T + efor£:>0, and 



thus, by Remark 15. 3[ boundedness and continuity of C on < t < T + e. □ 
Lemma 8.2. For M > and Ci> C » M sufficiently large, for 

Eo:= U1 + \x\Y'^\uo-u)\\h^ 

sufficiently small, and \S^'^^\bi + M|(5*"^^| < 2CEq, there exist solutions 
(m", 5^^, 5,") of (Km-lKm for all t > 0, satisfying 

(8.5) < CiEo 
and 

(8.6) + + M|5,"| < 2CEo. 
Proof. Define ( as in (18.41) . Then, it is sufficient to show that 

(8.7) at)<cEo+aiEo+my 

for fixed C, > 0, so long as the solution {u'\6'') G L^{t,H^{x)) x C°(t) 
of (l6ll]) -( l6:T2ll exists and 

(8.8) at) < (3/2)CEo, 
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in order to conclude that solution (m", 5") exists and satisfies fl8.8l) for all 
t >0, provided 

Eo < 25/2CC, 

is sufficiently small. 

For, by (KW and §1^, and ([23D, 

llw"(■,o)||^,5 = ||M^i^5 

II ~ —S II 
= \\Uo - U * Wffi 



n-1 



(8-9) < ||mo — w||//5 + —uWffi 

< Eo(l + 2ciC/M) 

is small, for Eq sufficiently small. Letting T be the maximum time up to 
which a solution (u", 5") exists and C ^ Co sufficiently small (note: T > by 
the weighted version of (18.91) . together with fl6.12p ). we find by Lemma [8?T| 
therefore, and the assumed bounds on and , that (m*^, 5") exists up 
to T + £, e > 0, and that ^ remains bounded and continuous up to T + e 
as well. Observing that (18.71) together with Eq < 2/(9C^ + 6C) implies that 
C(t) < (3/2) CEo whenever < (3/2) CEq, we find by continuity that 
C(t) < (3/2)^^0 up to t = T + e as claimed. 
By the definition of we obtain therefore 

(8.10) \u^\b, + |5"(t)(l + t)Uoo < {3/2)CEo. 

Thus, it is sufficient to establish first (18.71) . then afterward, assuming 

(8.11) \5^{t)\ < (CEo/4)(l + t)-i/2 
and 



(8.12) < {CEo/4M), 

from which we obtain (18.61) by summation with (18.101) : noting that (18.111) 
and (18.121) include the information that the integral equations for 5" and 5*" 
converge, we obtain also, by Lemma 16.11 and the fact that (m", 5") satisfies 
f l6ni|) - fl6:T2|) for all t > 0, that (m", 5", 5,") satisfies ([6S])-([6l]) as claimed. 
Finally, recalling (16. lip and applying Lemma l5.ll with 7 := 5" — 5"^^, we 
obtain (18. 5p so long as (18.71) remains valid, controlling Hu^H/^s by integrating 
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the righthand side of fl5.3l) and using fl8.8l) . the definition of (, and the as- 
sumed bounds on 5""^. (We carry out this last calculation in detail in the 
following paragraph, in the course of proving (18. 7p ). 

We now establish (18.71) assuming (I8.8p . By Lemma I5.H and the one- 
dimensional Sobolev bound |ti"|vi/3,oo < c\u"'\h5, we have 

\u-{t)\lr.<c\u^O)\U-'' 

+ c [ e'^'^^^-^^dM^lia + + |5"-5"-^n(r)dr 
Jo 



<c|M"(0)|^5e~^* 



(8.13) 



+ c y^e-'^^^^-^^dw^lia +max{|5"|M(5"-T})(r) dr 

< C2(^o'(l + ^ciC/Mf + (3C^o/2)') (1 + t)-'/' 

< (ClEo)'(l + ^)-'/^ 

for Ci > sufficiently large Eq sufficiently small, by (18.91) . (18.81) . and the 
definition of (. This verifies (18.51) . assuming (18. 8p . 

With (16. 2p . (18. 5p and the resulting Sobolev estimate llw^Hiys.oo < cCiEq, 
assumption \6\bi < '2CEq, and the definitions of ( and | ■ \bi, we obtain 
readily 



(8.14) 



\Q^'+R^' \ <c(C' + 4CX)(^ + '^'i)> 



and 

(8.15) l^^'l, \Sl*\<c{C + AC^El)^2, 

where and T defined in Lemmas I4.1H4.5[ 
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Expressing using fl6.6p as 



<(x,t) = / {H:'' + G:''){x,t;y)u^,-\y)dy 

J —oo 

Jo J -oo 

X S^'"'\6''~\6''~^){y,s)dyds 

"t p + oo 

H--\x, t - s; y)Sf 5'^-^){y, s)dyds 



J -oo 
t-1 



Gy^\x,t- s;y) 



—oo 

n-1 



(8.16) X iQ'*"~ +^'*"' iS^-\u-))iy,s)dyds 

nt r+oo 
J t—l J —oo 

X {Q'*''-\u-,u:)+R'*'"\5--\u^))^iy,s)dyds 

nt p+OO 

+ / / iH:-'-H;-'){x,t-s;y) 

Jo J -oo 

P+OO 

- / / H--\x,t-s;y) 

Jo J -OO 

X {Q'^''-\u-,ul) + R'''"\5^-\u-))Jy,s)dyds, 

and applying Lemmas 14.11 - 14.51 to (16.61) . (18.161) . and (16.121) . we thus obtain 
(18.71) as claimed. 

Likewise, we obtain easily (18. lip from (16. 8p and (18.90 . using Lemmas 
I4.1H4.5I and the definitions of C and | ■ \b^- 

Thus, it remains only to establish (I8.12p . This is more delicate, due to 
the appearance of M in the denominator of the righthand side, and depends 
on the key fact that estimate 5^^ of the asymptotic shock location is to 
linear order insensitive to the initial guess To see this, decompose the 
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expression (16.91) for 5*" into linear and nonlinear parts 

r+oo 



oo 

oo 



(8.17) 



J ~oo 

e|5.=o(y, +oo)(mo - u){y)dy 

e"-^ - e|5,=o)(2/, +oo)u'^~\y)dy 



oo 

oo 



oo 



Ia + h + Ic 



and 



11:=- / e'^-\y,+oo)S'^''-\S--\S^-')iy,s)dyds 

Jo J~oo 
p+oo n+oo 

- / / e"-i(y,+oo) 

Jo J -oo 

X {Q''''-\u\ul) + R'*^~\5--\u-))^{y,s)dyds, 

respectively. 

By estimates like the previous ones, we readily obtain 

\II\ < 2c(2CEo)^ 

which is << CEq/AM for Eq sufficiently small. Likewise, |/c| < cI^^I-Eq, by 
dHSD, flXT^ . and the Mean Value Theorem, hence is « CEq/AM for Eq 
sufficiently small (recall that we assume \5^\ < 2CEq), and 

\h\ < c\\uo - mIIli 

(8.19) <C2Ul + \x\')-'/\uo-mH^ 

< C2E0, 

hence is << CEq/AM for C > sufficiently large. 

Finally, Taylor expanding, and recalling fl2.4p and 03.131) . we obtain 

/+00 
e\s,=,{y,+^){du''' Id6,)\s,=,{y)dy 

18-^U) 



0{% 



I ii 
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which is also << CEq/AM for Eq sufficiently small (recall that we assume 
< 2CEq). Summing, we obtain fl8.12l) for Eq sufficiently small and C > 
sufficiently large, as claimed. This completes the proof. □ 



Proof of Theorem \1.4\ Define now 

(8.21) := \t\B, + M\j\, it,j) eB.xR. 



By Lemma [8.21 for M > sufficiently large, and 

and r > sufficiently small, T = {Ts, %J) is a well-defined mapping from 

S(0,r) C fii X M ^ fii X M. 

To establish the theorem, therefore, it suffices to establish that T is a contrac- 
tion on 5(0, r) in the norm | ■ |,,. For, then, applying Contraction Mapping 
Theorem, we find that 

(T5,T5j(5,5,) = (5,5,) 

has a unique solution (5", 5*") G 5(0, r) C -Bi x M, for which the associated 
(n", 5", 5,=") by Lemma [6.11 satisfy = u — u^*'^+^"i'^) with u a solution of 
(11.21) with initial data mq, and the stated decay estimates follow by (18. 5p and 

(ESD. 

That T is a contraction follows, provided we can establish on -6(0, r) the 
Lipshitz bounds 

(8.22) \r{5,5,) - r{lS.)l < a|((5,5,) - 
for some a < 1. 

Letting (u",5",5,") satisfy for and (u",5",5^) 

satisfy (I6.6I) - (I6.9I) with 5/"'^ replaced by 5"'"^, define variations 

(8.23) Au":=u"-u", A5" := 5" - 5", A^," := 5^ - 5^ 
and 

(8.24) A5"-i := 5"-^ - A(5r^ := - K'^- 
Likewise, define AG"~^, AH'"'^^, Ae""^ in the obvious way. 



29 



Differential variational equation. From fl6.1ip . we find after a brief calcu- 
lation that Am" defined by 



Au- = Am" - (m^"-^+^"-^W - u^-'+s-'i^)) 
satisfies the variational equations associated with 

(8.25) + f[u-)^ - {B{u-)u:)^ = (5"(t) - 6-\t)) — \,.-., 

from which we obtain by an energy estimate similar to that of Lemma 15.11 
and the observation 

the bound 

|AM"(t)|^4 < Ce-^*|AM"(0)|^4 

+ C f e-''^'-^\\AuX^ + max{|A(j"|2, |A(j"-i|2} 

(8.26) + I A^r^p max{|5"|2, |5"-^p})(r) dr 

< C Te"'^ (*-")( I Am" + |A(5"n(r)rfr 
+ C|(A5"-\A5ri)|^(l + t)-V2^ 

provided r (and so sup [[^"lliys and sup ||M"||iy5) is sufficiently small, so long 
as ||AM"||jy4 remains sufficiently small. We omit the (standard) details. 
Integral variational equations. Applying the quadratic Leibnitz formula 



30 



^ifg) = f^g + ^fg, we obtain 

/+0O 
{AH^-' + AG"-^)(x, t; y)ur\y)dy 
■oo 

/ + 00 
{H^-' + G^-'){x,t;y)Aur\y)dy 
-oo 



oo 

+ (Ai/"-i + AG'"-^)(x,t-s;y) 

Jo J -oo 

X S^*^~\5''-^,t-^){y,s)dyds 

t f+OO 

n— 1 I /^n—l\ 



+ / / {H--'-rG"'-'){x,t-s;y) 

Jo J -oo 

X AS{y, s)dyds 
(8.27) -/ / ^G;-\x,t-s;y) 



'0 -'-oo 

X {Q'*''-\u^,u2) + R'*'"\5^-\u^)){y,s)dyds 

ft f + OO 

G'^-'ix^t-s^y) 




J-oo 



X {AQ + AR){y,s)dyds 

rt r+oo 

+ / / AH--\x,t-s;y) 

Jo J-oo 

X {Q'*''-\u\u:) + R'*''-\S--\u^))^{y,s)dyds 
+ / / H--\x,t-s;y) 

Jo J-oo 

X {AQ + AR))^{y,s)dyds, 
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/oo 
-oo 

/oo 
er\y,t)i^ur\y)dy 
-oo 

ft z'+oo 

Aer'(y, i - s)S'*''-\5^-\ S^-'){y, s) dyds 

ft f + OO 

(8.28) ~ Jo J ^t~\y:t-s)AS{y,s)dyds 




'0 J~oo 

ft f + OO 



nt r+oo 

n-l 



^0 ^-oo 

X {AQ + AR{y,s)dyds, 

/ + 00 
Ae»-i(y,+ooX-i(y)ciy 
-oo 

/ + 00 
e^-\y,+oo)Au^-\y)dy 
-oo 

+ / / Ae''-\y,+^)S'^'''\6--\6''-'){y,s)dyds 

Jo J-oo 

p+oo r+oo 

(8.29) + / / e''-\y,+oo)AS{y,s)dyds 

Jo J-oc 

/• + 0O P + OO 

- / / Ae;-\y,+oo) 

Jo ./-oo 
Jo J-oo 



X 



(Ag + Ai?) {y,s)dyds, 
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and similarly for A 5", where 



(8.30) = (1^) 1..=.:'-^'^*"' + 0{\AS:-'\'e-^\-\ 

= OdA^r'e-^"!"!!). 

Now define 
(8.31) 

e(t):= sup f(|AM"| + |A<|)(^ + z^i+7A2)"'(a;,s) + |A5"(s)(l + s)|). 

x,0<s<t V / 

Let r' := |(A5"-i, A(5r^)|* = |A5"-i|si +M|A5:-^I, ^' sufficiently small. 
From fl8.26p and smallness of Iw'^liys and \u"'\h5, and the fact that r,r' « 1, 
we obtain 

(8.32) |Aw'^(-, t)|H4 < C{r' + e(t))(l + t)-^ 

which gives us a bound on L°°-norm of Am"^^, providing us, therefore, with 
the bounds, 

\AQ + AR\ < C{r^{t) + rr')(* + $i), 

(8.33) |Ag^ + ARy\, \AQyy + ARyyl < C {v ^ (t) + rr')T 

\AS\, \ASy\<Cir^{t) + rr')^2, 

Also, (187141) and (18151) hold with c(C^ + 4(:72e2) replaced with Cr^. We use 
(imj) . (I3J[5D and (I3J6D to obtain 

(8.34) Ae"-^ ~ eA(5, < r'e, 

and similar appropriate bounds for Aif""^, AG"~^ and their derivatives (of 
course, e in (I8.34p is defined at a point between and 5"~^). Next, using 
lemmas [4. 1H4.5I in a procedure parallel to the one used in the proof of Lemma 
18.21 we obtain 

m<c{r'+rm 

from which we conclude that 

Cr' 



1 - Cr 
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with constant C independent of r and r' . Now, replacing ^ in fl8.33l) with 
this bound, we plug back the result into (18.281) and into the similar formula 
for AJ". Notice that, with the exception of the first two terms, the other 
term in (18.281) have quadratic terms in their source term, so giving us small 
enough bounds. Hence, using one again lemmas I4.1H4.5[ we obtain 

(8.35) |A(5"i < {CE^r' + -^r' + Crr'){l + t)-\ 

of which the two first terms in the right hand side come from the first two 
terms of (I8.28p . Similarly we obtain 

(8.36) lA^'^l < (CEo + -^ + Cr)r'(l + t)-i 

We notice that {CEq + ^ + Cr) can be made arbitrarily small, provided that 
£^0) are small enough and M is large enough. Next, we use (I8.29P to bound 
using basically the same method used in (I8.17I) - (I8.20I) . and therefore 
obtaining 

(8.37) M\A6:\ < {CEo + Cry. 

This, together with (I8.35P and (18.360 . gives us (I8.22p with a < 1, finishing 
the proof of the (main) Theorem II. 4[ 

□ 

Remark 8.3. In order to control the norm of the variational problem, 
as in (I8.26p . we indeed need regularity C^ for the coefficients in hypothesis 
(HO), since one derivative is lost in variational energy estimate (15. 3p . 

Remark 8.4. HERE 

References 

[AMPZl] A. Azevedo, D. Marchesin, B. Plohr and K. Zumbrun, Nonuniqueness of 
Riemann solutions, Z. Angew. Math. Phys. 47 (1996), no. 6, 977-998. 

[Brl] L. Q. Brin, Numerical testing of the stability of viscous shock waves, Ph.D. 

dissertation, Indiana University, May 1998. 

[Br2] L. Q. Brin, Numerical testing of the stability of viscous shock waves, Math. 

Comp. 70 (2001) 235, 1071-1088. 



34 



[Br3] L. Brin, Numerical testing of the stability of viscous shock waves, Doctoral 

thesis, Indiana University (1998). 

[BrZ] L. Brin and K. Zumbnm, Analytically varying eigenvectors and the stability 

of viscous shock waves. Seventh Workshop on Partial Differential Equations, 
Part I (Rio de Janeiro, 2001). Mat. Contemp. 22 (2002), 19-32. 

[BDG] T. Bridges, G. Derks, and G. Gottwald, Stability and instability of solitary 

waves of the fifth-order KdV equation: a numerical framework, Phys. D 172 
(2002), no. 1-4, 190-216. 

[Pre] H. Frcistiihlcr, Some results on the stability of non-classical shock waves, J. 

Partial Diff. Eqs. 11 (1998), 23-38. 

[FreS] H. Frcistiihlcr and P. Szmolyan, Spectral stability of small shock waves. Arch. 

Ration. Mech. Anal. 164 (2002) 287-309. 

[PreZ] H. Frcistiihlcr and K. Zumbrun, Exam,ples of unstable viscous shock wa,ves, 

unpublished note, Institut fiir Mathcmatik, RWTH Aachen, February 1998. 

[GZ] R. Gardner and K. Zumbrun, The Gap Lemma and geometric criteria for 

instability of viscous shock profiles, Comm. Pure Appl. Math. 51 (1998), no. 
7, 797-855. 

[HZ] P. Howard and K. Zumbrun, Stability of undercompressive viscous shock 

waves, in press, J. Differential Equations 225 (2006), no. 1, 308-360. 

[HR] P. Howard and M. Raoofi, Pointwise asymptotic behavior of perturbed viscous 

shock profiles. Adv. Differential Equations 11 (2006), no. 9, 1031-1080. 

[HRZ] P. Howard, M. Raoofi, and K. Zumbrun, Sharp pointwise bounds for perturbed 

shock waves, J. Hyperbolic Differ. Equ. 3 (2006), no. 2, 297-374. 

[HuZ] J. Humphcrys and K. Zumbrun, Spectral stability of sm,all am,plitude shock 

profiles for dissipative symmetric hyperbolic-parabolic systems, Z. Angew. 
Math. Phys. 53 (2002) 20-34. 

[Kaw] S. Kawashima, Systems of a hyperbolic-parabolic composite type, with appli- 

cations to the equations of magnetohydrodynamics, thesis, Kyoto University 
(1983). 

[LZ2] T.P. Liu and K. Zumbrun, On nonlinear stability of general undercompressive 

viscous shock waves. Comm. Math. Phys. 174 (1995), no. 2, 319-345. 

[LRTZ] G. Lyng, M. Raoofi, B. Texier, and K. Zumbrun, Pointwise Green Function 
Bounds and stability of combustion waves, to appear, J. Diff. Eq. (2007). 



35 



[Ml] A. Majda, The stability of multi- dimensional shock fronts ~ a new problem 

for linear hyperbolic equations, Mem. Amer. Math. Soc. 275 (1983). 

[M2] A. Majda, The existence of multi- dimensional shock fronts, Mem. Amer. 

Math. Soc. 281 (1983). 

[M3] A. Majda, Compressible fluid flow and systems of conservation laws in several 

space variables, Sprmger-Verlag, New York (1984), viii+ 159 pp. 

[MP] A. Majda and R. Pego, Stable viscosity matrices for systems of conservation 

laws, J. Diff. Eqs. 56 (1985) 229-262. 

[MZl] C. Mascia and K. Zumbrun, Pointwise Green's function bounds and stability 

of relaxation shocks, Indiana Univ. Math. J. 51 (2002), no. 4, 773-904. 

[MZ2] C. Mascia and K. Zumbrun, Stability of small- amplitude shock profiles of 

symmetric hyperbolic-parabolic systems. Comm. Pure Appl. Math. 57 (2004), 
no. 7, 841-876. 

[MZ3] C. Mascia and K. Zumbrun, Pointwise Green's function bounds for shock 

profiles with degenerate viscosity. Arch. Rational Mech. Anal. 169 (2003), no. 
3, 177-263. 

[MZ4] C. Mascia and K. Zumbrun, Stability of large- amplitude shock profiles of 

hyperbolic-parabolic systems, Arch. Rational Mech. Anal. 172 (2004), no. 1, 
93-131. 

[MZ5] C. Mascia and K. Zumbrun, Stability of large- amplitude shock profiles of 

general relaxation systems, SIAM J. Math. Anal. 37 (2005), no. 3, 889-913. 

[PZ] R. Plaza and K. Zumbrun, An Evans function approach to spectral stability 

of small- amplitude viscous shock profiles, J. Disc, and Cont. Dyn. Sys. 10. 
(2004), 885-924. 

[Ra] M.-R. Raoofi, -asymptotic behavior of perturbed viscous shock profiles, J. 

Hyperbolic Differ. Equ. 2 (2005), no. 3, 595-644. 

[TZl] B. Texier and K. Zumbrun, Galloping instability of viscous shock waves, 



Preprint (2006); available http: //arxiv.org/abs/math. AP /0609331 



[TZ2] B. Texier and K. Zumbrun, Hopf bifurcation of viscous shock waves in 

compressible gas- and magnetohydrodynamics, preprint (2006); available 



http: / /arxiv.org/abs/math. AP /0612044( 



[ZH] K. Zumbrun and P. Howard, Pointwise semigroup methods and stability of 

viscous shock waves, Indiana Mathematics Journal V47 (1998) no. 4, 741-871. 



36 



K. Zumbruii. Multidirnensional stahiliiy of planar viscous shock waves, Ad- 
vances in the theory of shock waves, 307-516, Progr. NonUnear Differential 
Equations AppL, 47, Birkhuser Boston, Boston, MA, 2001. 

K. Zumbrun, Stability of large-amplitude shock waves for compressible 
Navier-Stokes equations, to appear, Handbook of Fluid Mechanics, Elsevier 
(2004). 

K. Zumbrun, Planar stability criteria for viscous shock waves of systems with 
real viscosity, to appear, CIME lecture notes series; preprint (2004). 

K. Zumbrun and D. Serre, Viscous and inviscid stability of multidimensional 
planar shock fronts, Indiana Univ. Math. J. 48 (1999) 937-992. 



37 



